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Abstract 

We discuss recently discovered links of the statistical models of normal random 
matrices to some important physical problems of pattern formation and to the 
quantum Hall effect. Specifically, the large N limit of the normal matrix model 
with a general statistical weight describes dynamics of the interface between two 
incompressible fluids with different viscousities in a thin plane cell (the Saffman- 
Taylor problem). The latter appears to be mathematically equivalent to the growth 
of semiclassical 2D electronic droplets in a strong uniform magnetic field with 
localized magnetic impurities (fluxes), as the number of electrons increases. The 
equivalence is most easily seen by relating the both problems to the matrix model. 
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1 Introduction 



The subject of the theory of random matrices is a random matrix M distributed with some 
probability measure d/i(M). Typically, one is interested in the distribution of eigenvalues 
and correlations between them as size of the matrix, N, tends to infinity. 

The range of physical applications of this theory is enormously vast, with the role 
played by the matrix M being very different. In complex systems or systems with disorder 
M turns out to be a good substitute for the Hamiltonian or transfer matrix. Physical 
characteristics of the system are obtained via averaging over one or another ensemble of 
large matrices. Most extensively employed and well-studied are ensembles of hermitian 
matrices. In more recent applications to statistical models on random lattices and to 
string theory the key element is a set of graphs in the diagrammatic expansion of random 
matrix integrals while M has no physical meaning by itself. For different aspects of 
random matrix theory and related topics see e.g. 0, [|, §]. 

Complex non-hermitian random matrices are employed in physics too. (A list of the 
relevant physical problems and corresponding references can be found in, e.g., [Ij.) New 
applications we are going to discuss are related to the distribution of their eigenvalues. 
To be specific, we consider the model of normal random matrices, i.e., such that M com- 
mutes with its hermitian conjugate, though similar results may hold for other ensembles. 
Eigenvalues of normal matrices are in general complex numbers. When N becomes large, 
they densely fill a domain in the complex plane, the support of eigenvalues, with the mean 
density outside it being zero. The shape of this domain is determined by the probability 
measure and by the size of the matrix. As N increases, the domain grows (see Fig. [I]). 
The growth law is our main concern in this paper. 

For simplicity we assume that the support of eigenvalues, D, is a connected domain. 
Let the size of the matrix grow linearly in time t: N ~ t. Then the support of eigenvalues 
grows in such a way that the normal velocity of the boundary is 

v(z) = grad</?(z) , z E dD (1) 

where ip is a function such that 

' Aif(z) = z e C\D 

< </?(z) ~ log |z| z — > oo (2) 

^ (p{z) = z G dD 

Here A = d 2 x + dy = 4:d z d 2 is the Laplace operator. We employ the complex notation 
z = x + iy. So, the dynamics of the boundary is governed by the function ip which is 
harmonic in the exterior of D with a source at infinity, and vanishes on the boundary. 
The solution of the boundary problem (0) is unique: <p(z) = \og\w(z)\, where w(z) is 
the conformal map from the exterior of the domain D onto the exterior of the unit circle 
such that oo is mapped to oo. Such a map exists by virtue of the Riemann mapping 
theorem. 

This dynamics is realized in various physical problems, classical and quantum. First, 
it is pertinent to viscous flows and other pattern formation processes when the normal 
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Figure 1: Growth of the support of eigenvalues. 



velocity of the moving front is proportional to the gradient of a harmonic field (see, 
e.g., ||). This mechanism is known as Darcy's law. The droplet of eigenvalues behaves 
like an incompressible fluid with negligible viscosity (say water) surrounded by a viscous 
fluid (say oil), two fluids being confined in a thin plane gap (the Hele-Shaw cell). Oil 
is withdrawn at infinity at a constant rate while water is injected. In this context, the 
function ip is identified with the pressure P in the viscous fluid with the opposite sign: 
ip = —P. In water, the pressure can be set to 0. The condition that P = on the 
interface amounts to neglecting the surface tension effects. This idealization is good 
until the curvature of the interface becomes large. When the surface tension is small 
enough, the dynamics becomes unstable. The moving interface develops many fingers, 
they split into new ones, and in a sufficiently long time the water droplet looks like a 
fractal. In the literature, this phenomenon is refered to as the Saffman- Taylor fingering. 
This growth process was linked to the matrix model in ||. 

As was recently pointed out in J7|, the same growth law applies to semiclassical 
dynamics of an electronic droplet confined in a plane on the lowest Landay level of a 
strong magnetic field. This suggests applications to the Quantum Hall effect. It turns 
out that the shape of the electronic droplet is sensitive to magnetic fluxes localized well 
away from it. As one changes degeneracy of the level to increase the number of electrons 
in the droplet, its shape evolves in accordance with the Darcy law, thus showing up 
fingering instabilities. This phenomenon is purely quantum. Like the Aharonov-Bohm 
effect, it is caused by quantum interference. The function ip in this case has no obvious 
physical interpretation. The characteristic scale of this phenomenon is less than that of 
the Saffman- Taylor fingering by a factor of 10 9 . Remarkably, the matrix model provides 
a unified mathematical treatment of the both phenomena. 

From mathematical side, it is also worth noting that calculation of certain expectation 
values and correlation functions of normal random matrices provides a constructive proof 
of some important mathematical statements in the inverse potential problem and the 
Dirichlet boundary problem proved in a different way in M. 

At last, the normal matrix model is known to be integrable. Its partition function is a 
tau-function of an integrable hierarchy of partial differential equations. Although we do 
not discuss integrability matters in this paper, let us point out that the above physical 
problems thus possess a hidden integrable structure. 
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2 Partition function of normal random matrices 



A matrix M is called normal if [M, M^] = 0, so that M and can be simultaneously 
diagonalized. The model of normal random matrices was introduced in [[|. The integrable 
structure of the model is quite similar to that of the more widely known model of two 
hermitian random matrices [T(| but physical interpretation is very different. 



The partition function is 

dfi (M)e^ v ^ M ^ (3) 



normal 



where h is a parameter and V is a real-valued function. The measure dfio is induced from 
the standard flat metric on the space of all complex matrices. To introduce coordinates 
in the subspace of normal matrices, one makes use of the decomposition M = UZU^ of a 
normal matrix M, where U is a unitary matrix and Z = diag (z\, . . . , z^) is the diagonal 
matrix of eigenvalues of the M. The measure is then given by 



i=i 



where dfj,u(N) is the Haar measure on the unitary group IA(N), and Aat(z) = Yu>j{ z i ~ z j) 
is the Vandermonde determinant. The partition function is, therefore, written as the 
following integral over eigenvalues: 



1 r N 

^N = ^J\A N (z)f]l(eT^d%) (4) 



(For notational simplicity we shall write V(z) instead of V(z, z).) This quantity has two 
important interpretations. 



If N 

One of them is the Coulomb gas picture |TT[. Writing Z N = — / Y\ d 2 Zj where 

i=i 



£= X>gi^-%l +h- l Y,V{ Zi ) (5) 

2D Coulomb energy potential 

we see that is the partition function of the 2D Coulomb gas in the external potential. 
Another one is the Quantum Hall picture suggested in 0. 



Quantum Hall picture. Consider spin-^ electrons on the plane in a non-uniform 
magnetic field B. The Pauli hamiltonian is 

H = — (Y-iW - A) 2 - ha 3 B) 

where A is the vector potential. If the magnetic field is uniform, the spectrum consists 
of equidistant Landau levels, each level being highly degenerate. The lowest level is very 
special. Due to a hidden supersymmetry of the problem, it can be found exactly and 
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remains highly degenerate even for arbitrary non-uniform field B JIj. The energy of 
this level equals while the degeneracy equals the integer part of the total magnetic flux 
$ = / B(z)d 2 z in units of the flux quantum $o — 27rfr (we set e = c = 1). One-particle 
states on the lowest level can be found explicitly. In the gauge A x = \d y V , A y = —\d x V 
they are 

= P„(*) <*p (^pj 

where B(z) = —hAV(z). Here P n = z n + [terms of lower degree] are holomorphic poly- 
nomials of any degree which is less than the degeneracy of the level [|T^j . 

Neglecting interactions between electrons, the wave function of iV particles on the 
lowest level is the Jastrow determinant: ~ tt^j det , n (^ m ), and so 



|*tf| 2 = ^|A*(*)| 2 e*£»^ (6) 

coincides with the statistical weight of normal random matrices expressed through eigen- 
values. The partition function (^) is, in this context, the normalization factor of the 
iV-particle wave function: / |\l/Ar| 2 IL d 2 Zi — Zn- The mean density of electrons coincides 
with the expectation value of the density of eigenvalues in the matrix model: 

. N-l 

NZ N l J \^ N {z, Ci, - - - , 6v-i)| 2 II d2 Zi = (tr«f (2) (* ~ M; 

i=l 

Similarly, multiparticle correlation functions are identified with multipoint correlation 
functions of densities. 

All the above relations are exact at any finite N. As h becomes small and iV large, one 
approaches a semiclassical regime. However, the semiclassical properties of the system 
are quite unusual. On the one hand, the density distribution acquires a well-defined 
edge, and one can speak about a well localized electronic droplet which behaves like an 
incompressible fluid. On the other hand, in this specific semiclassical regime, quantum 
effects are by no means negligible. Quite the reverse, they become rather strong if not 
dominant. In fact there is no surprize here because the semiclassical limit we are speaking 
about is not the usual one which would require excitations of higher energy levels. In 
our "semiclassical" limit all particles occupy the lowest level, so the droplet as a whole 
remains a quantum object. Amusingly enough, it is this limit where one makes contact 
with the purely classical Saffman- Taylor fingering. In the next section, we analyse the 
corresponding large N limit of the matrix integral. 



3 The semiclassical (large N) limit 

The large N limit we are interested in is iV — > 0, H — > with foiV finite and fixed. The 
expansion in iV -1 is then the same as the expansion in %. 

To elaborate the limit, we represent the energy £ (|5|) in the form 

£ = JJ p(z) log \z - z'\p(z')d 2 zd 2 z' + i J V{z)p{z)d 2 z 
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where p(z) = Y^i$ { z ~ z i) is the density of eigenvalues. In the limit, one treats p(z) as 
a continuous function normalized as / p(z)d 2 z = N. As % — > 0, both terms in (|5|) are of 
order N 2 , and the saddle point method can be applied to perform the integral. 

The saddle point condition is SS/Sp(z) = which yields the integral equation for the 
mean density: 

The meaning of this equation is especially clear in the Coulomb gas interpretation. It 
states that each charge is in the equilibrium. Indeed, consider a charge at the point z. 
The first term in the equation is the Coulomb force caused by other charges in the gas 
while the second term is the external force. The equation just tells that they compensate 
each other. Clearly, it makes sense to impose the equilibrium condition only in the 
domain where the charges are actually present, i.e., in the support of eigenvalues. 

So, the equation should be satisfied in a domain D where p ^ 0. Here we assume 
that D is a connected domain. For example, in the potential V = —\z\ 2 the eigenvalues 
uniformly fill the disk of radius VhN. Small perturbations of the potential slightly 
disturb the circular shape. 

It appears that in case of normal matrices the above integral equation is much easier 
to solve than the similar equation for distribution of eigenvalues of hermitian matrices. 
Indeed, on applying d z to both sides of eq. (^) we obtain p(z) = —j^rAV(z) in D, and 
p(z) = in C \ D. The domain D itself is determined by the condition 

d z ,V(z')dz' 

— — — — = 8 

dD Z — z' 

which can be derived from (|7|) with the help of the Cauchy integral formula. The condition 
means, in other words, that D is such that the function d z V on its boundary is the 
boundary value of an analytic function in C\D. 

An important particular case is V equal to — \z\ 2 plus a harmonic function which we 
parametrize by its Taylor coefficients tk at the origin: V = — \z\ 2 + 2TZe Y^tkZ k - Then 
the density is constant (equal to —A inside D and zero outside. The area of D is equal 
to UN . The shape of D is determined by the conditions 

1 f 

/ Z - k d 2 z = t k 

7TK JC\D 

which easily follow from (|8|), so that irktk are harmonic moments of the domain comple- 
mentary to D. In case of quadratic potential the domain is an ellipse [|l^ . 



The integrated version of eq. (0) allows one to find the leading contribution to the 
free energy, which is given by the value of h 2 £ at the saddle point. Let us denote the 
latter by F , then 



F » = -T^/i> Al/(2)loE 



1 1 

z z' 



AV(z')d 2 zd 2 z' (9) 



This is the tau-function of curves introduced in p. The leading asymptotics of the 
partition function as % — > is therefore — e F °^ h . 
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Figure 2: Normal displaycement of the contour. 



Small variations of the potential and the growth law. If one varies the potential, 
V — > V + SV, and size of the matrix, iV — > iV + 5N, the support of eigenvalues slightly 
changes its shape and area. Let us examine how it goes. It is natural to characterize an 
infinitesimal change of the boundary by its normal displaycement Sn(z) at each point z, 
so that Sn(z) is a continuous function on the boundary curve (see Fig. |[). 

First we vary the potential at constant N. The shape of the support of eigenvalues 
is determined by eq. @. Its variation can be written as 



z — z' 2 J z — z 

It is natural to employ the ansatz AV(z)Sn(z) = d n h(z) where h is yet unknown func- 
tion in the exterior of D such that h = on the boundary, and d n means its normal 
derivative, with the normal vector pointing outward. This ansatz is suggested by an 
easy transform of the second integral into a Cauchy integral. Combining the two terms, 
we get § 9 z'( 5V ( z J+ h ( z )) d z > = q f or z £ d This implies that the function d z (5V + h) is 
analytic in the exterior of D, i.e., that the function 5V + h is harmonic in there. Hence 
h(z) = 5V H (z) — 5V where, given a function /, we use the notation f H for the harmonic 
continuation of this function from the boundary to the exterior of D. Therefore, 



d n (5V H (z) - 5V(z) 

Similarly, to find how the domain D grows at constant V as N increases, we use 
eq. (|D again. This time the first term in ( |T0D is zero. Noting that 

5N = - — { 6n(z)£V(z)\dz\ 
Ann J 

it is easy to see that the solution is 

5n(z) = -^rd n log\w(z)\ (12) 

where w(z) is the conformal map from C \ D onto the exterior of the unit circle such 
that oo is sent to oo. Since log \ w\ = on the boundary, the normal derivative can be 
substituted by gradient. Therefore, at AV = const we get the Darcy law fll]) (cf. 0]). 
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a) b) 

Figure 3: a) The electronic droplet in the uniform magnetic held; b) The electronic 
droplet in the uniform magnetic held in the peresence of point-like fluxes at the points 
marked by <g>. 



Semiclassical electronic droplet in the presence of magnetic impurities. Let 

us apply the above results to the semiclassical behaviour of an electronic droplet in a 



strong magnetic field. The notion of the Quantum Hall droplet |14j implies that the 
electronic liquid is incompressible, i.e., all states at the lowest energy level are occupied. 
Therefore, we want the degeneracy of the level to be equal to N. This can be achieved in 
different ways. One of them is to assume the following arrangement. Let a strong uniform 
magnetic field B > be applied in a large disk of radius R . The disk is surrounded 
by a large annulus R < \z\ < R\ with a magnetic field B\ < such that the total 
magnetic flux through the system is NQ . The magnetic field outside the largest disk 
\z\ < R\ vanishes. The disk is connected through a tunnel barier to a large capacitor 
that maintains a small positive chemical potential slightly above the zero energy. If 
Bq is strong enough, the gap is large, and the higher levels can be neglected. In this 
arrangement, the circular droplet of iV electrons is trapped at the center. Its radius is 
much less than R . The function V(z) for \z\ < Ro is V(z) = —\Bq\z\ 2 . 

Now let us apply a non- uniform magnetic field SB somewhere inside the disk \z\ < Rq 
but well away from the droplet. Suppose that the nonuniform magnetic field does not 
change the total flux: / 5Bd 2 z = 0. The potential V(z) inside and around the droplet is 
modified as 

V(z) = - — \z\ 2 - - [\og\z- z'\5B(z')d 2 z' (13) 

2 7T J 

The second term is harmonic inside and around the droplet. One may have in mind 
thin solenoids carrying magnetic flux ("magnetic impurities"). In the case of point-like 
magnetic fluxes g; at points a iy we have V(z) = — ^B \z\ 2 + J2i Qi log \ z ~ a i\- 

In the presence of the fluxes, the shape of the droplet is no longer circular (Fig. ^) 
although the magnetic field inside the droplet and not far from it remains uniform and 
is not changed at all. In this respect this phenomenon is similar to the Aharonov-Bohm 
effect. The responce of the droplet to an infinitesimal change of the magnetic field 5B is 
described by eq. (ITT]) in which 



5V H (z) - 8V(z) = - [ G(z, z')5B(z')d 

71 JC\D 



'^ 2 z' 



3V 

Here G(z, z') is the Green function of the Dirichlet boundary problem in C\D normalized 
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in such a way that G(z, z') — > log \z— z'\ as z — > z'. In fact this formula holds for arbitrary 
55, not necessarily vanishing inside the droplet. In particular, for small point-like fluxes 
Sqi at some points a, t we have SV = Y,i ^Qi l°g | — «^ j , SB = —n Y,i SqiS^ 2 \z — ai), and 

SV H (z) - SV(z) = -J2 G(z, a t )S gi 

i 

If cii is inside, G(z, a^) is set to be zero. The sum, therefore, goes over outside fluxes only. 
The fluxes inside the droplet, if any, appear to be completely screened and do not have 
any influence on its shape. 

When B\ adiabatically increases, with Bo and SB fixed, the droplet grows because 
the degeneracy of the lowest level is enlarged and new electrons enter the system. The 
growth is described by eq. ( |T2"D with AV(z) = —2B which is equivalent to the Darcy 
law. 



4 Conclusion 

We have analysed the large N limit of the model of normal random matrices. It has been 
argued that as N increases, the growth of the support of complex eigenvalues simulates 
important physical phenomena: 

• Interface dynamics in viscous flows (the Saffman- Taylor fingering) in the zero sur- 
face tension limit 

• Semiclassical behaviour of 2D electronic droplets in the Quantum Hall regime 

The former is purely classical while the latter is purely quantum. 

The relation to the matrix model may help to suggest a way to regularize singular- 
ities which usually occur in the zero surface tension limit and to obtain an analytically 
tractable formulation of the Saffman- Taylor problem with surface tension. 
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